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ABSTRACT 


The spike (S) protein of coronaviruses, a type I membrane glycoprotein, is primarily responsible for entry into 
susceptible cells by binding with specific receptors on cells and mediating subsequent virus-cell fusion. The 
bovine coronavirus (BCoV) S protein is cleaved into two subunits, the N-terminal S1 and the C-terminal S2. 
The proteolytic cleavage site of S protein is highly conserved among BCoV strains and is located between 
amino acids 763 and 768 (KRRSRR). This study describes a single mutation in the S protein cleavage site of 
three Brazilian strains of BCoV detected in diarrheic fecal samples from calves naturally infected. The sequenced 
PCR products revealed that amino acid sequence of the cleavage site of our strains was KRRSSR, indicating 
a mutation at amino acid position 767 (R— S). This amino acid substitution occurred due to a single nucleotide 
substitution in the sequence of DNA corresponding to the proteolytic cleavage site, CGT to AGT. This is the 
first description of this nucleotide mutation (C to A), which resulted in the substitution of arginine to serine 
in the S cleavage site. In this study we speculated the probable effects of this mutation in the proteolytic 


cleavage site using the murine hepatitis coronavirus (MHV) as a comparative model. 
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Bovine coronavirus (BCoV), a member of the family 
Coronaviridae, order Nidovirales, belongs to group 2 of the 
coronaviruses which include murine hepatitis coronaviruses 
(MHV), porcine hemagglutinating encephalomyelitis virus 
(HEV), equine coronavirus (ECoV), rat coronavirus (RtCoV) and 
human respiratory coronavirus (HCoV-OC43) (9). BCoV is an 
enveloped virus with single-stranded, positive-sense RNA 
genome of approximately 32 kb length that encodes five major 
structural proteins: the nucleocapsid (N), the transmembrane 
(M), the hemaglutinin esterase (HE), the spike (S), and the small 
protein (E) (13). The S glycoprotein of BCoV is a large membrane 
glycoprotein of approximately 150 kDa that forms the peplomers 
(club-shaped structures) on the virion surface. The S protein is 
primarily responsible for the entry of coronavirus into 
susceptible cells by binding to specific receptors on cells and 
mediating virus-cell fusion and subsequent cell-cell fusion 
during infection (5). 


In several coronaviruses, such as infectious bronchitis virus 
(IBV), MHV and BCoV, as a late event in maturation, the protein 
is cleaved into two subunits: S1 (aminoterminal region) and S2 
(carboxyterminal region) (17). Proteolytic cleavage of the S 
protein of these coronaviruses occurs adjacent to a sequence 
of basic amino acids on the carboxyterminal region of S1. In the 
S protein of BCoV, a predicted basic amino acid sequence 
(KRRSRR) is involved in the cleavage by the host cell-derived 
proteolytic enzyme. This sequence, highly conserved among 
BCoV strains, encompasses amino acids 763 to 768; the cleavage 
occurs between amino acids 768 and 769 (19). The cleavage of 
the S protein has been reported as a process related to the viral 
infectivity and cell fusion from other group 2 coronaviruses. 
Studies related to MHV, the best-studied member of the 
Coronavirus family, demonstrated that cleavage of S is not 
essential for infectivity but is associated with enhanced cell 
fusion (syncytia) in infected cell monolayers (7,8,20). 
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This study describes a single mutation in the S protein 
cleavage site of wild type Brazilian strains of BCoV detected in 
calves naturally infected and speculates the possible biological 
effects of this mutation at the proteolytic cleavage site using 
the MHV as a comparative model. 

Three BCoV positive fecal samples (BR-UEL1, BR-UEL2 and 
BR-UEL3) were obtained from calves up to 30 days old with 
clinical signs of diarrhea in a Brazilian dairy herd from Minas 
Gerais State (21°41' 49" S; 45°18'45” W). These samples were 
previously identified as BCoV positive by RT-PCR assay for N 
gene detection (21) and negative for bovine group A rotavirus 
and Cryptosporidium sp by polyacrylamide gel electrophoresis 
technique and modified Ziehl-Nielsen method. Fecal samples 
were diluted 2-fold in 0.01 M phosphate-buffered saline (PBS) 
pH 7.2 (137 mM NaCl; 3 mM KCl; 8 mM Na,HPO,; 15 mM 
KH.HPO,), centrifuged at 3000 x g for 15 min at 4°C to remove the 
cell debris. The supernatants were used for RNA extraction. 
Aliquots of 400 pl from fecal suspensions were treated with SDS 
at a final concentration of 1% (w/v), homogenized by vortexing 
and kept at 56°C for 30 min. A combination of phenol/chloroform/ 
isoamyl alcohol and silica/guanidinium isothiocyanate methods 
was performed according to Barreiros et al. (1) with slight 
modifications. Briefly, 400 ul of phenol/chloroform/isoamyl 
alcohol (25:24:1) were added, vortexed and heated at 56°C for 15 
min (16). The mixture was centrifuged at 10,000 x g for 10 min and 
the supernatant was transferred into a new tube and processed 
by the silica/guanidinium isothiocyanate method (2). The RNA 
was eluted from the silica pellet with 50 ul of diethyl- 
pyrocarbonate (DEPC) treated sterile water by 15 min incubation 
at 56°C and centrifuged at 10,000 x g for 10 min. The supernatant 
fraction was kept at -20°C until further use. 

Specific oligonucleotide primers that flank the sequence 
corresponding to S protein cleavage site were designed using 
the Gene Runner program version 3.05 (Hastings Software Inc., 
Hastings, NY) (http://www.generunner.com). The primers 
sequences (positions calculated from the start codon of the S 
gene) were: SPK7_F: 5’-TAACTCTTCCGAACCAGCA-3’ (nt 
2085-2103) and SPK7_R: 5’-AATCGCTTCCTAAACAACC-3’ 
(nt 2701-2719), which amplify a predicted fragment of 636 bp. 

The reverse transcription (RT) reaction was performed with 
a mixture of 8 pl of extracted RNA and 2 ul of the random primer 
pdN6 (GE Healthcare, Little Chalfont, UK) and incubated at 
97°C for 4 min. Subsequently, it was placed on ice for 5 min and 
10 pl of RT mix containing 1x RT buffer (50 mM Tris-HCl, pH 8.3, 
3 mM MgCl, 75 mM KCl), 0.1 mM of each dNTP (Invitrogen™ 
Life Technologies, USA), 10 mM DTT, 100 units of M-MLV 
Reverse Transcriptase (Invitrogen™ Life Technologies, USA) 
and ultrapure sterile water to a final volume of 20 ul were added 
and incubated at 37°C for 60 min and followed by inactivation at 
95°C for 5 min. 

For amplification, 8 pl of the RT reaction were added to 42 ul 
of the PCR mix containing 1.5 x PCR buffer (30 mM Tris-HCl pH 
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8.4and75mM KCl), 2mM MgCh, 0.2 mM of each dNTP, | ul (20 
pmol) of each primer, 2.5 units Platinum Taq DNA polymerase 
(Invitrogen™ Life Technologies, USA) and ultrapure sterile 
water to a final volume of 50 ul. The reaction was performed in 
a thermocycler (PTC-200, MJ Research Co. Water Town, MA, 
USA) with the following time and temperature conditions: one 
step of 4 min/94°C; followed by 40 cycles of 1 min/94°C, Imin/ 
52°C, 1 min/72°C and a final step of 7 min/72°C. 

The products were analyzed by electrophoresis in a 2% 
agarose gel in TBE buffer pH 8.4 (89 mM Tris; 89 mM boric acid; 
2 mM EDTA), stained with ethidium bromide (0.5 ug/ml) and 
visualized under UV light. 

The PCR amplicons were purified using GFX PCR DNA and 
Gel Band Purification (GE Healthcare, Little Chalfont, UK) and 
sequenced in MegaBACE 1000/Automated 96 Capillary DNA 
Sequencer (GE Healthcare), according to the manufacturer’s 
instructions. Sequencing was performed in both directions using 
the forward (F) and reverse (R) primers. Sequences quality 
analysis was performed using Phred/Phrap/Consed Analysis 
Program (http://www.phrap.org) and sequence similarity 
search was performed using BLAST software (http:// 
www.ncbi.nlm.nih.gov/BLAST/). The nucleotide sequences of 
the wild-type Brazilian BCoV strains were aligned and compared 
with the BCoV reference strains using the CLUSTAL/W program. 

The sequences reported in this study have been deposited 
in the GenBank database under accession numbers: 
DQ479421 (BR-UEL1), DQ479422 (BR-UEL2) and DQ479423 
(BR-UEL3). The predicted secondary structure, hydrophilicity 
profile and protease map were determined using Protean of 
DNAStar software (DNAStar, Inc., Madison, WI, USA). 

The amino acid sequence of the cleavage site of the three 
wild-type Brazilian BCoV strains was KRRSSR, indicating an 
amino acid mutation at position 767 (R— S). Divergence in amino 
acid sequence of the S protein cleavage site among 
coronaviruses strains are indicated in Fig. 1. 

The only mutation described in the proteolytic cleavage site 
of bovine coronaviruses was reported in 1991 with the wild type 
French BCoV strain F-15, which also showed an amino acid 
substitution at position 767 (RV) (Fig. 1). However, the authors 
did not discuss the effects related to the amino acid change (24). 

Sequence analysis suggested that the BCoV BR-UEL 
sequences were trustworthy and not the result of RT-PCR or 
sequencing errors, since the same mutation was successfully 
amplified and sequenced in all samples and also in different 
fecal aliquots from each animal. Furthermore, the mutations 
associated with viral adaptation in cell culture were excluded, 
since the samples did not undergo cell culture passage prior to 
PCR amplification and sequencing. 

Although there are no study related to mutational analysis 
in the proteolytic cleavage site of BCoV strains, this 
phenomenon has been previously described for MHV strain. 
Computational analysis demonstrated that the HCoV-OC43, 
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Figure 1. Predicted amino acid sequences for the spike gene in the cleavage site proximities. Residues in bold represents proteolytic 
cleavage signal sequence of the Mebus strain. The arrow between amino acids 768 and 769 indicates the cleavage site between S1 
and S2. Shaded residue indicates the mutation in the cleavage signal. BCoV BR-UEL strains (GenBank accession numbers 
DQ479421, DQ479422, DQ479423), Mebus (U00735), BCoV ENT (AF391541), LSU94 (AF058943), OK5 14 (AF058944), Kakegawa 
(DQ479424), LY 138 (AF058942), Quebec (AF220295), F15 (D00731), Vaccine (M64668), L9 (M6466). 


MHV-AS59 and BCoV presented high degree of identity among 
group 2 coronaviruses (18). Based on these evidences, the MHV 
was used as an appropriate model to hypothesize the effects of 
the new mutation observed in the S protein cleavage signal of 
wild-type BCoV strains during this study. 

Cleavage of S protein is related to efficient cell-cell fusion 
(syncytia) by MHV in infected cell monolayers but is not 
necessary for virus-cell fusion (infectivity). The MHV spike 
proteins, that have mutations that eliminate cleavage into S1 
and S2 subunits, carry out cell-to-cell fusion very inefficiently; 
however, they are able to mediate entry into susceptible cells 
with similar efficiency as wild-type virus (3,8, 14). Although the 
cleavage of the S protein is not a prerequisite for fusion, it does 
enhance the induction of syncytia (3). Hingley et al. (12) and 
Hann et al. (10) reported that the cleavage of the S protein of 
MHV-AS59 mutant was greatly reduced in comparison with the 
wild-type virus. The sequence at the predicted cleavage site of 
MHV-AS59 wild-type is RRAHR while that of MHV-A59 mutant 
is RRADR. The MHV-A59 mutant showed a substitution of a 
weakly basic histidine (H) by negatively charged aspartic acid 
(D) in the cleavage signal (8). These authors suggest that the 
introduction of a negatively charged amino acid into this highly 
basic region destroys the signal or otherwise prevents its 
recognition and thereby inhibits cleavage of the spike. 


Similar to MHV, BCoV S gene also has a proteolytic cleavage 
site formed by a group of basic amino acids upstream from the 
S2 amino terminal region. Except for BCoV F15, this sequence 
is highly conserved in all known BCoV strains (KRRSRR). In 
the BCoV BR-UEL samples was identified an amino acid 
substitution (RS) at amino acid position 767 (KRRSSR) (Fig. 
1). This amino acid change occurred due to a single nucleotide 
mutation at the proteolytic cleavage site, CGT to AGT. This 
mutation (C to A), that resulted in the change of the arginine 
codon to serine, had never been described. The cleavage into 
the two subunits occurs between amino acid 768 and 769 and 
is thought to be mediated by cellular trypsin-like proteases 
(19). The protease map indicated that three BR-UEL strains 
lack one predicted cleavage site by trypsin when compared 
with other known BCoV. In MHV studies there are strong 
evidence that furin enzyme is the cell protease responsible for 
cleavage of the S protein in cultured cells (10). In addition, 
although it did not alter the predicted secondary structure or 
hydrophilicity of the protein, the substitution of arginine to 
serine is a nonconservative substitution. Nonconservative 
amino acid changes within a single protein may result in 
alterations of the physical energy of the protein-protein 
interaction and may destabilize the protein native conformation 
(6). Navas-Martin et al. (15) suggest that the proper spatial 
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arrangement of the S1 and S2 subunits is crucial for the 
biological functions of the S protein. 

Yoo and Deregt (23) generated BCoV mutants and confirmed 
that a single point mutation at domain II of $1 protein was 
responsible for the escape of BCoV from immunological 
response. Similarly, He et al. (11) described that a single amino 
acid substitution in the receptor-binding domain of SARS 
coronavirus S protein disrupted the antigenic structure and 
binding activity. In addition, a single amino acid change has 
already been demonstrated to influence MHV ability to spread 
within the central nervous system (22). 

Recently, Navas-Martin et al. (15) have demonstrated that a 
single cleavage signal substitution in MHV mutants may play a 
major role in the virulence. These authors have associated this 
mutation with higher viral load and highly virulent phenotypes. 
Interestingly, the BR-UEL strains were obtained during an 
outbreak of neonatal diarrhea that culminated with the death of 
calves in a Brazilian dairy cattle herd. Because BCoV infections 
resulting in high mortality is uncommon, the results described 
by Navas-Martin et al. (15) for MHV strain should be also 
investigated for BCoV proteolytic cleavage site mutants. 
Divergence within the S1 gene of others BCoV Brazilian strains 
had already been described by Brandao et al. (4). These authors 
reported the first description of a gap of 18 nucleotides (deletion 
of 6 amino acids) within the hypervariable region of the S1 
subunit from their isolates and that was also found in human 
coronavirus strain OC43. 

We are uncertain if this single mutation prevented the normal 
proteolytic cleavage of the S protein or increased viral load, as 
occurred in MHV mutants. Therefore, additional experiments 
must be done to determine the effects of the amino acid change 
(RS) within the proteolytic cleavage site on the structural 
and functional features of the BCoV S protein. Development of 
a system which introduces infectious cDNA clones with specific 
mutations into the BCoV genome will provide an important tool 
to determine the role of these BCoV mutants in the pathogenesis 
of the neonatal calf diarrhea. 
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RESUMO 


Identificacéo de uma mutacao no sitio de clivagem da 
proteina da espicula em amostras brasileiras de 
coronavirus bovino 


A proteina da espicula (S), uma glicoproteina de membrana 
do tipo I, é primariamente responsavel pela entrada do virus em 
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células susceptiveis por meio da interacgdo inicial com receptores 
celulares especificos e subseqiiente mediacao da fusao virus- 
célula. A proteina S do coronavirus bovino (BCoV) é clivada em 
duas subunidades: a S1, na regiao N-terminal e a S2, na regiao 
C-terminal. O sitio de clivagem proteolitica da proteina S é 
altamente conservado entre as estirpes de BCoV e esta situado 
entre os aminoacidos 763-768 (KRRSRR). Este estudo descreve 
uma mutacao no sitio de clivagem da proteina S de trés estirpes 
do BCoV detectadas em amostras fecais diarréicas de bezerros 
naturalmente infectados no Brasil. O seqiienciamento dos 
produtos de PCR identificou a seqiiéncia de aminodcidos 
KRRSSR no sitio de clivagem de nossas amostras, indicando 
uma muta¢do na posi¢4o 767 (RS). Esta mutacdo ocorreu 
devido a uma Utnica substituigéo de nucleotideo no sitio de 
clivagem proteolitica, alterando o cédon CGT para AGT. Esta é 
a primeira descri¢ao desta muta¢ao de nucleotideo (C para A), 
que resultou na substituigaéo do aminodcido arginina por serina 
no sitio de clivagem da proteina S. Neste estudo também sao 
sugeridos os provaveis efeitos desta mutacdo no sitio de 
clivagem proteolitica utilizando o coronavirus da hepatite dos 
camundongos (MHV) como um modelo comparativo. 


Palavras-chave: BCoV, seqiienciamento, proteina da espicula, 
sitio de clivagem, mutagao 


REFERENCES 


1. Barreiros, M.A.; Alfieri, A.F.; Médici, K.C.; Leite, J.P.; Alfieri, A.A. 
(2004). G and P genotypes of group A rotavirus from diarrhoeic 
calves born to cows vaccinated against the NCDV (P[1],G6) rotavirus 
strain. J. Vet. Med. B, 51, 104-109. 

2. Boom, R.; Sol, C.J.A.; Salimans, M.M.M.; Jansen, C.L.; Wertheim- 
van Dillen, P.M.E.; Van Der Noordaa, J. (1990). Rapid and simple 
method for purification of nucleic acids. J. Clin. Microbiol., 28, 
495-503. 

3. Bos, E.C.W.; Heijnen, L.; Luytjes, W.; Spaan, W.J.M. (1995). 
Mutational analysis of the murine coronavirus spike protein: effect 
on cell-to-cell fusion. Virology, 214, 453-463. 

4. Brandao, P.E.; Gregori, F.; Richtzenhain, L.J.; Rosales, C.A.R.; 
Villarreal, L.Y.B.; Jerez, J.A. (2006). Molecular analysis of Brazilian 
strains of bovine coronavirus (BCoV) reveals a deletion within the 
hypervariable region of the S1 subunit of the spike glycoprotein 
also found in human coronavirus OC43. Arch. Virol., 151, 1735- 
1748. 

5. Cavanagh, D. (1995). The coronavirus surface glycoprotein. In: 
Siddell, S.G. (ed.), The Coronaviridae, Plenum Press, New York, 
USA, p. 73-113. 

6. Davis, N.L.; Fuller, F.J.; Dougherty, W.G.; Olmsted, R.A.; Johnston, 
R.E. (1986). A single nucleotide change in the E2 glycoprotein gene 
of Sindbis virus affects penetration rate in cell culture and virulence 
in neonatal mice. Proc. Natl. Acad. Sci. USA, 83, 6771-6775. 

7.  Frana, M.F.; Behnke, J.N.; Sturman, S.; Holmes, K.V. (1985). 
Proteolytic cleavage of the E2 glycoprotein of murine coronavirus: 
host-dependent differences in proteolytic cleavage and cell fusion. 
J. Virol., 56, 912-920. 

8. Gombold, J.L.; Hingley, S.T.; Weiss, S.R. (1993). Fusion-defective 
mutants of mouse hepatitis virus A59 contains a mutation in the 
spike protein cleavage signal. J. Virol., 67, 4504-4512. 


10. 


11. 


12. 


13%, 


14. 


tS. 


16. 


Gonzales, J.M.; Gomez-Puertas, P.; Cavanagh, D.; Gorbalenya, A.E.; 
Enjuanes, L. (2003). A comparative sequence analysis to revise the 
current taxonomy of the family Coronaviridae. Arch. Virol., 148, 
2207-2235. 

Hann, C.A.M.; Stadler, K.; Godeke, GJ.; Bosch, B.J.; Rottier, P.J.M. 
(2004). Cleavage inhibition of the murine coronavirus spike protein 
by a furin-like enzyme affects cell-cell but not virus-cell fusion. J. 
Virol., 78, 6048-6054. 

He, Y.; Li, J.; Jiang, S. (2006). A single amino acid substitution 
(R441A) in the receptor-binding domain of SARS coronavirus spike 
protein disrupts the antigenic structure and binding activity. Biochem. 
Biophys. Res. Commun., 344, 106-113. 

Hingley, S.T.; Bombold, J.L.; Lavi, E.; Weiss, S.R. (1994). MHV- 
AS59 fusion mutants are attenuated and display altered hepatotropism. 
Virology, 200, 1-10. 

Lai, M.M.C.; Cavanagh, D. (1997). The molecular biology of 
coronaviruses. Adv. Virus. Res., 48, 1-100. 

Leparc-Goffart, I.; Hingley, S.T.; Chua, M.M.; Phillips, J.; Lavi, E.; 
Weiss, S.R. (1998). Target recombination within the spike gene of 
murine coronavirus mouse hepatitis virus-A59: Q159 is a determinant 
of hepatotropism. J. Virol., 72, 9628-9636. 

Navas-Martin, S.; Hingley, S.T.; Weiss, S.R. (2005). Murine 
coronavirus evolution in vivo: functional compensation of a 
detrimental amino acid substitution in the receptor binding domain 
of the spike glycoprotein. J. Virol., 79, 7629-7640. 

Sambrook, J.; Fritsch, E.F.; Maniatis, T. (1989). Molecular cloning: 
a laboratory manual, 2nd Ed. Cold Spring Harbor Laboratory Press, 
Cold Spring Harbor, New York. 


17. 


18. 


20. 


21. 


22. 


23. 


24. 


Spike protein cleavage site 


Spann, W.; Cavanagh, D.; Horzineck, M.C. (1988). Coronaviruses: 
structure and genome expression. J. Gen. Virol., 69, 2939-2952. 
St-Jean, J.R.; Jacomy, H.; Desforges, M.; Vabret, A.; Freymuth, F.; 
Talbot, P.J. (2004). Human respiratory coronavirus OC43: genetic 
stability and neuroinvasion. J. Virol., 78, 8824-8834. 


. Storz, J.; Rott, R.; Kaluza, G. (1981). Enhancement of plaque 


formation and cell fusion of an enteropathogenic coronavirus by 
trypsin treatment. Infect. Immun., 31, 1214-1222. 

Sturman, L.S.; Ricard, C.S.; Holmes, K.V. (1985). Proteolytic 
cleavage of the E2 glycoprotein of murine coronavirus: activation 
of cell-fusing activity of virions by trypsin and separation of two 
different 90K cleavage fragments. J. Virol., 56, 904-911. 

Takiuchi, E.; Stipp, D.T.; Alfieri, A-F.; Alfieri, A.A. (2006). Improved 
detection of bovine coronavirus N gene in faeces of calves infected 
naturally by a semi-nested PCR assay and an internal control. J. 
Virol. Methods, 131, 148-154. 

Tsai, J.C.; Groot, L.; Pinon, J.D.; Iacono, K.T.; Phillips, J.J.; Seo, 
S.H.; Lavi, E.; Weiss, S.R. (2003). Amino acid substitutions within 
the heptad repeat domain 1 of murine coronavirus spike protein 
restrict viral antigen spread in the central nervous system. Virology, 
312, 369-380. 

Yoo, D.; Deregt, D. (2001). A single amino acid change within 
antigenic domain II of the spike protein of bovine coronavirus 
confers resistence to virus neutralization. Clin. Diagn. Lab. Immunol., 
8, 297-302. 

Zhang, X.; Kousoulas, K.G.; Storz, J. (1991). Comparison of the 
nucleotide and deduced amino acid sequences of the S genes specified 
by virulent and avirulent strains of bovine coronaviruses. Virology, 
183, 397-404. 


703 


